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SUMMARY 


Rate  equations  for  the  densities  of  free  excitons  and 
excitons  bound  to  two  types  of  neutral  acceptors  in  silicon  are 
solved  for  steady  state  in  the  absence  of  saturation.  These  rate 
equations  explicitly  include  the  terms  for  reverse  and  forward 
tunneling  of  bound  excitons  from  one  type  of  neutral  impurity 
to  another.  The  tunneling  rates  are  calculated  using  a  simple 
model  of  an  exciton  in  a  one-dimensional  semi-infinite  double 
potential  well .  The  energy  eigenvalue  equation  for  an  exciton 
in  this  potential  well  is  derived  for  estimating  the  exciton 
tunneling  time.  The  steady-state  solutions  of  the  rate  equations 
yield  an  expression  for  the  ratio  of  the  bound  exciton  lumines¬ 
cence  intensity  as  a  function  of  the  impurity  concentrations. 

The  relative  photoluminescence  intensities  for  the  systems 
Si:(B,In),  Si:(Al,In),  Si:(Ga,In),  Si:(B,Al),  Si:(B,Ga),  and 
Si(Al,Ga)  are  calculated  for  various  values  of  the  relevant  free 
exciton  capture  cross  section  ratios.  This  model  predicts  no 

exciton  tunneling  for  any  of  the  above  systems  for  the  low 
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impurity  concentration  range  of  10  -10  cm  .  For  the  systems 

with  large  differences  in  the  bound  exciton  energy  levels  such  as 

Si:(B,In),  Si:(Al,In),  and  Si:(Ga,In),  and  having  indium  concen- 
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trations  exceeding  10  cm  ,  it  predicts  quenching  of  shallow 
impurity  bound  exciton  luminescence  because  the  forward  exciton 
tunneling  rate  from  the  shallow  level  to  the  deep  level  of  indium 
dominates  and  the  reverse  exciton  tunneling  rate  from  indium  to 
the  shallow  impurities  is  negligible.  For  the  systems  with  small 
differences  in  the  bound  exciton  energy  levels  such  as  Si:(B,Al) 
and  Si:(B,Ga),  the  theory  predicts  enhancement  of  shallow  impurity 
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bound  exciton  luminescence  beyond  certain  concentrations  depending 
upon  the  free  exciton  capture  cross  section  ratios  because  in 
these  cases  the  reverse  exciton  tunneling  rate  dominates.  For 
tne  system  Si:(Al,Ga)  in  which  the  difference  in  the  bound  exciton 
energy  levels  is  very  small,  gallium  bound  exciton  luminescence 
dominates  when  the  gallium  concentration  exceeds  lO^cm  ^  if  the 
aluminum  free  exciton  capture  cross  section  is  less  than  the 


FOREWORD 


Materials  requirements  for  silicon-based  devices  are 
becoming  ever  more  stringent  as  VHSIC-related  applications  become 
more  important  to  Air  Force  research.  Photoluminescence  (PL)  has 
been  shown  to  be  an  effective  method  for  impurity  and  defect  con¬ 
centration  determination  in  the  10^-10^cm  ^  range.  However, 
recent  experimental  and  theoretical  evidence  indicates  that  the 

PL  method  may  have  limited  usefulness  as  a  quantitative  technique 

15  -  3 

for  impurity  concentrations  above  10  cm  .  This  paper  estab¬ 
lishes  a  theoretical  framework  for  examining  the  exciton  tunneling 
problem  in  doubly-doped  silicon  and  thus  establishes  a  "window  of 
applicability"  for  quantitative  concentration  measurements  in 
silicon. 

One  of  the  authors  (D.  S.  Moroi)  would  like  to  thank  the 
AFOSR  for  supporting  the  preliminary  part  of  the  work  under  the 
1982  USAF-SCEEE  Summer  Faculty  Research  Program,  Contract  No. 
F49620-82-C-0035  and  Universal  Energy  Systems  for  funding  the 
project  for  the  summer  months  of  1983  under  Contract  No. 
F33615-82-C-5001  and  Task  No.  729-074.  He  would  like  to  extend 
his  appreciation  to  the  Laser  and  Optical  Materials  Branch, 
Electromagnetic  Materials  Division,  the  Materials  Laboratory  of 
the  AFWAL  for  providing  him  with  the  opportunity  to  spend  the 
very  worthwhile  and  rewarding  summers.  It  is  our  pleasure  to 
thank  Dr.  Patrick  M.  Hemenger  for  providing  useful  references  and 
inspiring  discussions.  His  collaboration,  guidance,  and  con¬ 
tinuing  encouragement  during  the  course  of  research  is  appre¬ 
ciated.  We  are  also  thankful  to  Dr.  Frank  Szmulowicz  for 


numerous  stimulating  discussions  on  modelling  exciton  tunneling. 
Finally,  we  want  to  thank  Mr.  David  O'Quinn  for  his  fine  work  in 
producing  many  of  the  figures  for  the  tunneling  factors  and 
Ms.  Niki  Maxwell  for  a  fine  job  of  typing  the  manuscript. 
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SECTION  I 
INTRODUCTION 


Characterization  of  semiconductor  materials  in  terms  of 
their  electromagnetic  properties  is  important  to  the  development 
of  semiconductor  devices  and  various  electrical  and  optical  tech¬ 
niques1  have  been  utilized  to  study  those  electromagnetic  proper¬ 
ties.  Particularly,  recent  progress  in  the  technology  of  inte¬ 
grated  circuits  demands  near-perfect  crystals  of  silicon  and 
requires,  therefore,  a  great  improvement  in  the  methods  of 
analysis  for  minute  amounts  of  impurities  and  defects  in  silicon 
since  they  drastically  change  its  electromagnetic  properties. 

Recently,  attempts  have  been  made  to  use  photoluminescence 

(PL)  as  a  quantitative  tool  for  measuring  impurity  concentrations 
,  2 

Tajima  obtained  boron  and  phosphorus  concentrations  in  silicon 

crystals  by  PL  analysis  at  liquid-helium  temperature.  Nakayama 

et  al . 3  measured  and  analyzed  the  low  temperature  exciton 

luminescence  of  silicon  on  the  basis  of  rate  equations  describing 

the  formation  and  decay  kinetics  of  free  excitons  (FE) ,  bound 

excitons  (BE)  and  bound  multiexciton  complexes.  Tajima  and 
4 

Nomura  applied  the  PL  technique  to  the  quantitative  analysis 
of  shallow  impurities  incorporated  intentionally  and  unintention¬ 
ally  in  silicon  epitaxial  layers.  Mitchard  and  McGill5  used  PL 
to  determine  relative  concentrations  in  the  Si:(B,In)  system 
while  Brown  et  al.  investigated  the  optical  and  electrical  prop¬ 
erties  of  the  systems  Si:(In,Al)  using  PL,  Hall  effect  transport, 
infrared  absorption,  and  photoconductivity  measurements.  In  both 
papers,  it  was  observed  that  the  luminescence  from  the  shallower 
acceptor  (B  or  Al)  of  a  two-acceptor  system  may  be  quenched  when 


the  concentration  of  the  deeper  acceptor  (In)  exceeds  certain 

critical  values  (>10^cm  ^).5-6  In  other  words,  no  aluminum 

nor  boron  bound  exciton  luminescence  was  observed  despite  high 

concentrations  of  A1  and  B  in  the  samples  used.  It  is  believed 

that  the  forward  tunneling  of  bound  excitons  to  In  bound  states 

is  responsible  for  the  quenching  of  the  B  and  A1  luminescence. 

Consequently,  it  appears  that  PL  may  not  be  particularly  useful 

for  determining  the  concentration  of  a  second  impurity  in  the 

presence  of  another  dopant  at  much  higher  concentrations. 

The  main  goal  of  this  research  project  is  to  understand 

the  electromagnetic  properties  of  multiply-doped  silicon,  and,  in 

particular,  to  provide  a  theory  interpreting  the  experimental 

£ 

results  on  PL  from  the  Si:(In,Al)  system  by  Brown,  et  al .  ,  as 
well  as  predicting  the  reverse  exciton  tunneling  for  the  systems 
Si:(B,Al),  Si:(B,Ga),  and  Si:(Al,Ga).  We  do  not  attempt  to  solve 
the  PL  problem  using  the  exact  Schrodinger  equation  for  an  exciton 
because  of  its  complexity;  instead,  we  use  the  simplest  model  to 
obtain  the  PL  intensity  ratio  for  the  doubly-doped  semiconductor, 
leaving  more  complicated  and  sophisticated  models  for  later 
investigations . 

In  the  next  section,  we  present  a  general  PL  theory  for 
doubly-doped  semiconductors;  this  theory  takes  into  account  for¬ 
ward  and  reverse  exciton  tunneling,  impurity  distributions,  the 
resulting  rate  equations,  and  their  steady-state  solutions  in 
the  absence  of  saturation.  We  calculate  in  Section  III  the  for¬ 
ward  and  reverse  exciton  tunneling  rates  using  a  simple  model  of 

a  BE  in  a  one-dimensional  semi-infinite  double  potential  well 

7 

( OSDW )  using  standard  methods  of  quantum  mechanics  .  The  exciton 


2 


energy  eigenvalue  equation  is  also  derived  for  estimating  the 
tunneling  time.  In  Section  IV,  we  calculate  the  tunneling  rate 
and  derive  the  energy  eigenvalue  equation  for  a  BE  in  a  one¬ 
dimensional  finite  double  potential  well  (OFDW)  for  future  use. 
The  numerical  results  and  discussion  along  with  recommendations 
and  suggestions  for  further  work  are  presented  in  Section  V. 


SECTION  II 


A  THEORY  OF  PHOTOLUMINESCENCE  FOR  A  DOUBLY-DOPED 
SEMICONDUCTOR  INCLUDING  EXCITON  TUNNELING 

1.  IMPURITY  DISTRIBUTION 

To  study  the  concentration  effect  on  luminescence  of  two 
impurity  atoms  in  a  semiconductor,  we  consider  here  some  direct 

g 

interactions  between  an  exciton  bound  to  impurity  atom  "1"  (A^) 
and  its  nearest  neighbor  neutral  impurity  atom  "2"  (A^)  which 
results  in  tunneling  of  the  exciton  from  the  atom  A^  to  the  atom 
(due  to  overlap  of  the  wavef unctions  for  the  A^  BE  and  BE) 
and  its  reverse  tunneling.  From  here  on,  we  will  use  the  sub¬ 
script  and  "2"  for  quantities  associated  with  impurity  atom 

"1"  with  shallow  BE  energy  levels  and  with  impurity  atom  "2" 
with  deep  BE  energy  levels,  respectively;  e.g.,  N^  =  the  concen¬ 
tration  of  B  atoms,  I2  =  the  intensity  of  PL  from  the  In  BE, 

=  the  FE  capture  cross  section  by  impurity  atom  "1". 

We  calculate  first  the  probability  that  the  nearest- 
neighbor  atom  A2  is  between  r  and  r  +  dr  from  a  given  atom  A^ 
as  shown  in  Fig.  1 

?2(r)dr  =  4Tir2drN2  exp  [-^irr3N2 1  (1) 

where 

2 

4nr  drN2  =  The  probability  that  an  atom  A2  is  between  r  and 
r  +  dr  from  a  given  atom  A^ 

4  3 

exp[-yT;r  N0]  =  The  probability  that  there  is  no  atom  A ^  between 
the  atom  A^  and  r  (obtained  from  the  Poisson 
distribution) . 
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A^  =  Impurity  Atom  "l1 


o-----r--->  I© 


No  A2  atoms 


A^X  =  Impurity  Atom  "2' 


Figure  1.  Configuration  of  impurity  atoms  in  a 
doubly-doped  semiconductor. 


Note  that 


P2(r)dr  =  1 


The  average  value  of  r  for  a  given  A^  concentration  is 


<r>2  = 


/  rP2(r)dr  =  (3/4ttN2)  1/3  F  ( 4/3) 


where  T(x)  is  the  gamma  function9.  For  N2=1016cm~3  we  find 

'r>2  =  257°.  Considering  that  the  Si  BE  Bohr  radius  is  43°,  it  is 

reasonable  to  assume  that  direct  interaction  effects  are  important 

Now  the  density  of  atoms  A^  with  a  nearest-neighbor  atom 
A2  between  r  and  r+dr  is  given  by 


Nx(r)  =  NjP2(r) 


so  that 


OO  Z*00 

J0  Nl(r)dr'Nl/o  P2(r)dr  *  N1 


Simply  switching  the  subscript  "2"  to  subscript  in  Eqs .  (1) 

and  (4),  we  obtain  the  probability  that  the  nearest-neighbor 
atom  A^  is  between  r  and  r+dr  from  a  given  atom  A 2  as  well  as 
the  density  of  A2  atoms. 
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2. 


THE  RATE  EQUATIONS  AND  STEADY  STATE  SOLUTIONS 


We  use  a  model  for  the  formation  and  decay  kinetics  of 
FE  and  BE  including  forward  and  reverse  BE  tunneling  to  derive 
a  formula  for  the  relative  intensity  of  PL  from  two  dopants  in 
a  semiconductor  as  a  function  of  their  concentrations.  For  this 
purpose,  we  start  with  the  rate  equations  for  the  densities  of 
the  FE  and  BE8. 


dnFE 

dt 

9  (VFE 

+  Y1  +  Y2)nFE  +  °lnl  + 

p2n2  ' 

(6) 

-Jn1  ( r) 

3 1 

Yl(r)nFE  - 

{ (r)  }n1  (r) 

+  0.2  (r)  n2 

(r)  , 

(7) 

>n2 (r) 

3 1 

Y2(r)nFE  * 

{  v2  +  p2  +  il2  (r)  }n2  (r) 

+  &1(r)n1 

(r)  . 

(8) 

Here 


n 


FE 


FE 


n  . 
D 


ni  (r) 


(r) 


(r) 


=  the  free  exciton  (FE)  generation  rate 
=  the  FE  density 

=  the  FE  decay  rate 

=  the  density  of  excitons  bound  to  impurity  atom 
j  (j  =  1  or  2);  e.g.,  "1"  =  A1 ,  ,,2,,  =  In 
=  the  density  of  excitons  bound  to  impurity  atom  i 
with  a  nearest-neighbor  impurity  atom  j  i) 
between  r  and  r  +  dr 

=  the  FE  capture  rate  by  impurity  atom  i 

=  the  FE  capture  rate  by  impurity  atom  i  with  nearest- 

neighbor  impurity  atom  j  (^  i)  between  r  and  r  +  dr 
=  the  BE^  decay  rate 

=  the  thermal  release  rate  of  BE. 

l 

=  the  distance  dependent  exciton  tunneling  rate  from 
impurity  atom  i  to  impurity  atom  j  i) 


Integrating  Eq.  (13)  over  r,  we  obtain 


r 

n .  =  /  n . 

1  Jn  1 


( r)  dr 


Y4 


=  q-FE 11  -  ^  jj  -  q1 


where 


J  . 


1  -  X 


-1 


Pj(r)w^(r)  {1  +  w^(r)  +  w2(r)}  dr 


with 


P.(r)  =  4tir^N,  exp  ( -^y-r  ^N_.  ) 


THE  INTENSITY  RATIO  OF  BOUND  EXCITON  LUMINESCENCE 
Using  Eq .  (17)  we  have  the  ratio  of  PL  intensity  of 


BE^  to  that  of  FE 


I  . 
l 


“FE 


f  .  n  .  f  .  v  . 

l  =  ri  c.u  +  12  J.  -  j.} 

fFLnFE  fEF  1  Yi  3  1 


where  f.  is  the  BE.  oscillator  strength  and  f_„  is  the  FE 
ii  3  FE 

oscillator  strength.  With  the  help  of  Eq.  (19) ,  we  have  the 
intensity  ratio  of  BE  luminescence 


N, 


=  R 


12  N. 


where 


R 


12 


R  ^2 }  Rd'2) 


wi  th 


R 


(o: 


12 


i  XX2  ( 
I  -..2N1  ' 


^  o1/(v1  +  p.^ 


f2  o2/  (\>2  +  P2) 


-  the  value  of  R^2  in  the  absence  of 


exciton  tunneling, 


R(l,2) 


(23) 


and 


1  +  n*J1  -  J2 


[1 

V2 


°1N1 

°2^2 


(24) 


We  name  R(l,2)  in  Eq .  (23)  the  "tunneling  factor"  of  a  bE  from 

impurity  1  to  impurity  2.  It  has  the  following  characteristics: 
no  BE  tunneling  will  occur  if  R(l,2)  =  1,  the  forward  BE 
tunneling  from  the  shallow  impurity  to  the  deep  impurity  dom¬ 
inates  if  R(l,2)  <  1,  and  the  reverse  BE  tunneling  from  the  deep 
impurity  to  the  shallow  impurity  dominates  if  R(l,2)  >1.  We  point 

out  here  that  Eq.  (21)  contains  the  vital  information  on  PL  from 
a  doubly-doped  semiconductor  which  depends  on  the  impurity  con¬ 
centrations  as  well  as  parameters  such  as  the  FE  and  BE  oscillator 
strengths  and  the  FE  capture  cross  sections. 


SECTION  III 


TUNNELING  RATE  AND  ENERGY  EIGENVALUE  EQUATION 
FOR  AN  EXCITON  IN  A  ONE  DIMENSIONAL 
SEMI-INFINITE  DOUBLE  POTENTIAL  WELL 

Before  we  compare  this  theory  with  the  experimental  results, 
we  must  calculate  the  exciton  tunneling  rates  12^  (r)  in  Eqs .  (16)  , 
(18),  and  (23).  The  general  topic  of  tunneling  in  solids  has  been 
considered  extensively  in  the  literature10,11. 

1.  TUNNELING  RATE  CALCULATION 

Here  we  treat  the  exciton,  a  correlated  electron-hole  pair, 
as  a  single  particle  and  assume  that  it  moves  in  a  OSDW,  as  shown 
in  Fig.  2.  It  is  assumed  that  the  exciton  is  initially  in  poten¬ 
tial  well  "1"  (Region  I  in  Fig.  2) . 


Figure  2.  An  exciton  in  one-dimensional  semi-infinite  double 
potential  well  (OSDW) . 


The  motion  of  an  exciton  with  effective  mass  m  in  the 


potential  well  shown  in  Fig.  2  is  governed  by  the  Schrodinger 
equation 


*  d  2  ^ 

_  _  +  {E  -  v),  =  o 


(25 


with  the  solutions 


Y z  ( x )  =  AelklX  +  Be  lklX 


0  <  x  <  a. 


,,,  ,  .  „  ~KX  ^  KX 

V 1 1  ( x )  =  Ce  +  De 


a^  <  x  <  a^  +  b 


(26 


¥HIU)  -  Ee 


ik2x 


a^  +  b  <  x  <  a^  +  a^  +  b 


where 


ki  =  /2m(V.  -  E)/fi 


=  /2mE/fi 


(27 


m  =  0.6m 


and  constants  A,  ...E  are  to  be  determined  by  the  boundary  con 
ditions.  The  boundary  conditions  between  regions  I  and  II, 
and  II  and  III  yield 

i  ( 


E/A  =  4  i  .'^e  ‘ 


/A. 


=  e  rb(l  +  ij(l  +  U2)  -  eKb(l  -  i^)  (1  -  ic2)  (28 


t  .  —  k  .  ti  . 

1  11 


<  =  k2  ( a^  +  b)  ,  c,^  =  k^/K 


A  transmission  coefficient  T  can  be  calculated  by  taking 
the  ratio  of  the  current  density  of  particles  transmitted 
through  the  barrier,  J  j. ,  and  the  incident  current  density,  J. 


11 


(29) 


T  =  J  /Jj  =  (k2/k1) (E/A) " 

=  rk^k2<2/A 
A  =  K- I A1 I 2/4 

=  (k^2  +  <2)(k22  +  k2 ) sinh2 (Kb)  +  <2(k^2  +  k22) 

The  tunneling  rate  from  well  "1"  to  well  "2”  is  defined  as 
the  transmission  coefficient  divided  by  the  period  of  classical 
oscillations  within  potential  well  "1": 


r  =  a1  +  a2  +  b 

In  arriving  at  Eq.  (30) ,  we  have  ignored  the  term  with  e 


in  Eq.  (29)  assuming  a  large  separation  of  two  wells.  The  reverse 


tunneling  rate  is  obtained  by  exchanging  subscripts  l«-*2  in  Eq.  (30) 


multiplied  by  a  Boltzmann  factor  exp  (-AE/kT)  with  AE  =  -  E2 

(the  difference  in  the  BE  energy  levels  in  two  wells. 


The  tunneling  rate  from  impurity  atom  A^  to  impurity  atom 

A0  by  Eq.  (30)  is  overestimated  for  the  approximate  concentration 

17  -3 

range  below  N2  =  10  cm  for  the  actual  three-dimensional  problem. 
We  discuss  below  how  to  remedy  this  problem.  It  is  reasonable  to 


assume  that  12  impurity  atoms  A2  are  distributed  at  equi-distance 
<c^2  from  each  other  on  the  spherical  shell  of  radii  <r>2  an<^ 
<r-_+dr  centered  at  the  atom  A.,  which  is  analogous 


to  a  close  hexagonal  packing  of  atoms  in  a  solid.  The  total 

width  of  the  potential  wells  of  the  12  impurity  atoms  on  the 

shell  does  not  occupy  the  entire  region  of  the  spherical  shell, 

if  the  concentration  of  the  atoms  is  less  than  approximately 
17  -3 

10  cm  .  Taking  account  of  this  fact,  the  more  accurate 


expression  for  the  tunneling  rate  is  given  by 


for  N_  <  1017cm'3  (32) 

Z  ft 

C  =  1  for  N„  >  1017cm"~3 

1  2  S3 


17  -3 

For  concentrations  greater  than  about  10  cm  ,  the  correction 
factor  becomes  unity  due  to  overlap  of  the  potential  wells. 

The  relationship  between  the  depth,  width,  and  the  energy 
eigenvalue  for  the  exciton  in  an  isolated  semi-infinite  well  is 
given  by 

tan(k.a.)  +  k./*-  =  0  (3  3) 

ii  l 

We  assume  that  the  reverse  tunneling  from  well  "2"  to 
well  "1"  will  take  place  after  the  exciton  in  well  "2"  is  raised 
from  its  deeper  energy  level  E^  to  a  shallower  energy  level 
by  thermal  excitation. 

The  tunneling  rate  may  be  calculated  by  various  other  means 
e.g.  the  semiclass ica 1  approximation  by  Wentzel ,  Kramers,  and 

7 

Brillouin  (the  WKb  approximation)  yields  a  tunneling  rate 
similar  to  Eq.  (10).  In  general,  the  tunneling  rate  depends 


on  a  dominant  factor  exp(-2<r)  regardless  of  the  model  used 
because  the  wavef unction  representing  a  BE  behaves  as  exp(-<r) 
after  tunneling  through  a  potential  barrier. 


2.  DERIVATION  FOR  THE  ENERGY  EIGENVALUE  EQUATION 

Here  we  derive  the  energy  eigenvalue  equation  for 
an  exciton  in  an  OSDW  in  order  to  compute  the  BE^  tunneling  time, 
The  wavef unctions  in  the  three  regions  which  satisfy  the  appro¬ 
priate  boundary  conditions  are 

l^I  =  A'sin(k1x) 

Vjj  =  B 1  e-KX  +  C'eKX  (34) 

Yj-jj  =  D'sin(k2(a1  +  a2  +  b  -  x)  ) 

By  eliminating  the  constants  A’,  B',  C',  and  D'  from  the 
four  equations  for  the  boundary  conditions  (4'1(a^)  =  fjjla^), 

Vj'  (ax)  =  YII'  (ax)  ;  V TJ  +  b)  =  V  ^  (a±  +  b)  ,  4^'  (ax  +  b)  = 
rIIl'^al  +  b)  ) »  we  have  the  exciton  energy  eigenvalue  equation 


F1F2 


where 


-2Kb 


tan (k . a . )  +  k . /< 
ii  r 

tan (k . a . )  -  k . /< 

ii  i' 


Equation  (35)  reduces  to  Eq.  (33)  for  a  large  separation,  b,  of 
the  two  wells. 
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SECTION  IV 


TUNNELING  RATE  AND  ENERGY  EIGENVALUE  EQUATION 
FOR  AN  EXCITON  IN  ONE  DIMENSIONAL  FINITE 
DOUBLE  POTENTIAL  WELL 


We  give  here  the  tunneling  rate  of  an  exciton  in  a  OFDW 
and  the  energy  eigenvalue  equation  for  future  use.  The  exciton 
tunneling  rate  is 


i  <r>  =  qe“2Kr 


‘1 

r 


=  4e 


-k (ai  +  a  2 ) 


(3‘ 


j(al  +  a2  ^  +  ^ 


Tne  realistic  tunneling  rate  is  obtained  by  multiplying  Eq.  (37) 
by  the  correction  factor  given  in  Eq.  (32)  as  in  the  case  for 
the  OSDW . 

The  energy  eigenvalue  equation  is 

-2kd 


G1G2 


=  e 


(38) 


where 


G. 


1 2 cot  (k.  a.  )  -  (—  -  A-)| 


/  h 


k  . 
i 


SECTION  V 

NUMERICAL  RESULTS  AND  DISCUSSIONS 

It  is  convenient  to  introduce  the  forward  tunneling 

coefficient  w.  (=  w. .)  which  is  defined  by 
1  i]  J 

w.  =  w..  =  Q . / (v .  +  p.)  (40) 

i  ij  ii  i 

where  (v.  +  p^)  is  the  total  decay  rate  of  BE^ ,  i.e.,  the  sum 
of  the  BE^  decay  rate  and  thermal  release  rate  of  BE^ ,  each  of 
which  is  estimated  using  the  experimental  value  of  the  BE  life- 

5  12 

time  at  4.2K.  '  It  turns  out  that  v. ,  the  thermal  release 

l 

rate,  is  negligible  for  all  cases  with  the  exception  of  the  BE 
for  boron  at  this  temperature.  We  have  used  the  experimental 
values  of  ionization  energies  determined  by  Lipari  et  al.13 
to  estimate  the  BE  ionization  energies  with  the  help  of  the 
Haynes  rule.  There  are  still  four  parameters  a^ ,  ,  and 

V2  (the  widths  and  depths  of  the  two  potential  wells)  to  be 
determined  before  calculating  the  tunneling  coefficients;  however 
the  energy  eigenvalue  equation  for  each  separate  well  (Eq.  (33)) 
effectively  reduces  the  number  of  unknowns  to  two.  For  the 
system  Si:(B,In),  we  have  fixed  the  numerical  value  of  one 
parameter,  a  ,  to  match  the  values  of  R(B,In)  obtained  by  Mitchard 

D 

and  McGill5  for  the  boron  concentration  range  Nn  =  10^-10^cm  3 

D 
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and  the  indium  concentration  range  NJn  =  10  -10  cm  and  approx 

imated  the  other,  aIn,  as  the  first  Bohr  radius  in  the  context 
of  the  hydrogenic  BE  model.  For  the  rest  of  the  systems,  we 
have  determined  the  shapes  of  the  potential  wells  using  a  similar 
process.  The  numerical  values  of  the  BE  ionization  energy  (E) , 
wavenumber  (k),  potential  depth  (V),  potential  width  (a),  total 
decay  rate  (v+p),  forward  tunneling  coefficient  (w. .)  and  reverse 


tunneling  coefficient  (w^)  for  each  system  are  summarized  in 
Table  1  for  the  BE  model  presented  in  Section  III  (an  exciton 
in  a  OSDW) . 

Note  here  that  the  reverse  tunneling  coefficients  w. 

3i 

are  calculated  with  the  help  of  Eq.  (40)  exchanging  the  sub¬ 
script  i  j  multiplied  by  the  Boltzmann  factors.  As  seen  in 
Table  1,  the  reverse  tunneling  coefficients  w_.^  for  the  systems 
Si:(B,In),  Si:(Al,In)  and  Si:(Ga,In)  are  extremely  small  and  thus 
neglected  for  evaluating  their  tunneling  factors  R(i,j). 

With  few  exceptions,  the  FE  capture  cross  sections  by 

neutral  acceptors  in  silicon  are  not  well  known.  We  take 
5  14 

o^/jin  =0.1  '  for  the  system  Si:(B,In).  For  the  system 

Si:(Al,In),  we  have  estimated  a  . /a  =  1.4,  using  a  linear 

a  i  in 

relationship  between  the  reciprocal  of  the  FE  capture  time15 
and  the  FE  capture  cross  section.  However,  we  also  examine  four 
otner  values  0.05,  0.2,  5,  and  20  for  o,../a_  to  evaluate  the 
tunneling  factor  R(Al,In)  because  of  uncertainty  in  the  FE 
capture  cross  section  ratio.  We  use  the  same  set  of  values  for 
the  o_  /o  rate  in  the  system  Si:(Ga,In)  as  that  for  the  system 
Si:(Al,In)  because  aluminum  and  gallium  have  similar  electronic 
properties  in  Si.  For  the  systems  Si:(B,Al)  and  Si:(B,Ga),  we 
take  the  set  of  values  0.005,  0.02,  1/14,  0.5,  2,  and  20  for 
their  FE  capture  cross  section  ratios.  Finally,  for  the  system 
Si:(Al,Ga),  we  examine  values  of  0.05,  0.2,  1,  2,  10,  and  50  for 


the  0^/0  ratio. 


ions  for  the  tunneling  factor 

are  plotted  in  Fig.  3a,  as 

of  Nb  (1012-1020  cm-3) . 

m t  views  of  tne  three  dimensional 

N„  and  NT  .  Similar  plots  for 

by  Mitchard  and  McGill5  are 

m.  For  the  boron  concentration 

the  indium  concentration  range 

ineling  factors  R(B,ln)  and 

each  other  because  the  potential 

Si  :  (B , In)  have  been  chosen  to 

i  above  concentration  ranges. 

ive  ranges,  however,  R(B,In) 

ises  more  rapidly  than  Rmm(B,In) 

:e  tunneling  coefficient  of 
7 

that  of  10  by  Mitchard  and 

in  the  exponential  factor  for  the 

In  other  words,  the  present 

ling  rate  and  thus  stronger 

luminescence  than  that  predicted 
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centrations  below  10  cm  ,  both 
ity  which  implies  no  BE  tunneling 

.  3b,  3c,  4b,  and  4c  indicate 
(B , In)  =  1)  . 

\)  for  the  system  Si:(Al,In) 


5a-9c.  For  all  the  values  of 

A] 

no  BE  tunneling  i f  both  aluminum 
14  -3 

less  than  10  cm  .  Comparison  c 
=  0.05)  and  9a  (R(Al,In) 
there  is  more  BE  tunneling  at  a  g 
larger  value  of  •  ,/u_  .  In  Tabl 
R  ^  (Al,In)  and  the  tunneling  fac 
'  . /o  for  N  .  =  3.03  x  lO^cm  3 
in  order  to  compare  the  theory  an< 
et  al.6  The  values  of  R^°^  (Al,I) 
using  the  relative  oscillator  str< 
et  al.  The  theory  predicts  vei 


TABLE 

R(0)(Al,In)  AND  R(Al,In) 
FOR  Na1  =  3.03xl015cm“3  Ms 


0  /  o 

AY  In 

R(o)  (Al , In ) 


0.05 

0.0875 


k  (Al ,  In)  0 . 000 :■! 7 


0.2 


0.350 


0 . UU0334 


factors  for  the  above  concentratio 
the  help  of  Eqs .  (20)  and  (21) ,  th 
the  system  Si:(Al,In)  is 


Referring  to 


relative  PL 


Table  2, 
intensity 


R (o) (Al , In 
Eq.  (41)  giv 


cross  section  ratios  of 
!0  and  are  plotted  in  Figs. 


20 


3.57  x  10 


(42) 


which  is  still  below  the  sensitivity  of  the  experiment  to  observe 
PL  from  aluminum.  Thus  it  confirms  the  experimental  results  by 
Brown  et  al.6  that  no  aluminum  BE  luminescence  should  be  observed 
at  these  concentrations. 

The  tunneling  factors  R(Ga,In)  for  the  system  Si:(Ga,In) 

are  calculated  for  the  o_  / oT  =  0.05,  0.2,  1.4,  5,  and  20  and 

Ljo.  jl  n 

plotted  in  Figs.  10a-14c.  These  figures  are  very  similar  to 

those  for  the  system  Si:(Al,In),  i.e.  they  indicate  no  BE 

tunneling  from  gallium  to  indium  if  the  concentrations  are  below 
14  -3 

10  cm  ,  and  strong  quenching  of  gallium  bound  exciton  lumines- 

17  -3 

cence  for  In  concentrations  beyond  10  cm 

For  the  systems  Si:(B,Al)  and  Si:(B,Ga),  we  have  calculated 

the  tunneling  factors  R(B,Al)  and  R(B,Ga)  for  the  FE  capture  cross 

section  ratios  °B/°Al  =  0B/°Ga  =  0.005,  0.02,  1/14,  0.5,  2,  and 

20,  and  plotted  them  in  Figs.  15a-20c  and  Figs.  21a-26c, 

respectively.  The  figures  corresponding  to  the  same  FE  capture 

cross  section  are  very  similar  to  each  other.  Once  again  they 

show  no  BE  tunneling  from  the  shallow  impurity  to  the  deep 

14  -  3 

impurity  if  their  concentrations  are  below  10  cm  .  For 

~  0.005,  however.  Figs.  15a  and  21a  show  very 
little  quenching  of  boron  bound  exciton  luminescence  and  rather 


14  -3 

an  enhancement  if  Nr,  '  10  cm  ,  N,. 

B  A1 


lO^cm  ^  and  NT 


t  i nlo  -3 
3x10  cm 


The  finite  values  of  tne  reverse  tunneliny  coefficients  are  respon¬ 
sible  for  this  phenomenon  (see  Table  1) .  In  other  words,  for  a 
small  value  of  the  FE  capture  cross-section  ratio,  the  reverse 
tunneling  from  the  deeper  impurity  to  the  shallower  impurity 
starts  to  dominate  the  forward  tunneling  from  the  shallower 


I 


21 


impurity  to  the  deeper  impurity  for  the  above  concentrations, 
resulting  in  only  slight  quenching  of  the  shallow  impurity  bound 


exciton  luminescence.  For  o  /u  =  o/o  -  20,  Figs.  20a  and 

B  A1  B  Gel 


26a  show  considerable  quenching  of  boron  bound  exciton  lumines¬ 


cence  if  N_.  <  lO^crvt  ^  and  N.,  =  N„  =  BxlO^^cm  ^  and  enhance- 
B  A1  Gel 


,14  -3 


merit  if  <  10  cm  and  the  extremely  high  concentrations  of 


20  -3 

=  N  =10  cm  .  For  a  large  value  of  the  FE  capture  cross¬ 


section  ratio,  a  large  number  of  excitons  are  captured  by  the 
shallower  impurity,  most  of  them  are  lost  to  the  deeper  impurity 
via  forward  BE  tunneling,  and  the  smaller  number  of  excitons 


captured  by  the  deeper  impurity  stay  there  because  of  the  small 
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reverse  BE  tunneling  coefficients,  if  N  <  10  cm  and  Na1 

B  AX 


16  —  3 

N„_  =  3x10  cm  .  Thus  considerable  quenching  of  boron  bound 


Ga 


exciton  luminescence  occurs  at  these  concentrations.  For  the 


20  -  3 

extremely  high  concentrations  of  =  NGa  =  10  cm  ,  the  reverse 


tunneling  from  the  deeper  impurity  to  the  shallower  impurity 
starts  to  dominate  the  forward  tunneling  and  enhances  the 
shallower  impurity  bound  exciton  luminescence.  Once  again,  the 
plateau  regions  of  Figs.  15  a, c  -  26b, c  indicate  the  regions 


exhibiting  no  BE  tunneling. 

The  results  of  the  calculations  for  the  tunneling  factors 
E(Al,Ga)  for  the  system  Si: (Al , Ga)  with  the  FE  capture  cross- 
section  ratios  ••  ./•  =  0.05,  0.2,  1,  2,  10,  and  50  are  plotted 

A  i.  Gel 

iii  Figs.  27a  -  32c.  Figures  27a  and  28a  show  no  quenching  of 
aluminum  bound  exciton  luminescence  and  its  enhancement  if 
.  ,  •  lO^cm  ^  and  N„  >  3xl0^cm  ^  duo  to  the  fact  that  the 

t  V.i  C j 

lorwnrci  and  reverse  tunneling  coefficients  are  of  the  same  order 
of  magnitude  and  the  difference  in  their  BE  energy  levels  is  very 


small.  For  o  /.;  ■  1,  slight  quenching  of  aluminum  bound 

AX  CjcI 

exciton  luminescence  occurs  at  certain  gallium  concentrations; 

e.g.  N  =  3xl0^^cm  ^  for  a  ../o_  =  1  and  N„  =  3xlU^cn\  ^ 

Oci  A  L  Gel  (ju 


f°r  'o,/ 


Al'  ~  Ga 


50  (see  Figs.  29a  and  32a) .  Figs.  29a  and  32a 


also  show  enhancement  of  aluminum  bound  exciton  luminescence 
if  N'  ■  lO^^cm  ^  for  oA1/o  =  1  and  N„  >  6xl0^^cm  ^  for 
j,,/^  =  50,  because  the  reverse  tunneling  starts  to  dominate 
beyond  these  concentrations. 

Finally,  to  ensure  the  occurrence  of  BE  tunneling,  we  have 
estimated  the  tunneling  time  of  a  BE  from  one  potential  well  to 
another  for  the  system  Si:(B,In)  at  NT  =  lO^cm  ^  by  two  differ- 
ent  methods: 

1.  The  reciprocal  of  the  tunneling  rate  (Eq.  (30)) 
gives  the  tunneling  time  of  1.14x10  ^ s . 

2.  Eq.  (35)  with  the  help  of  time-dependent  per¬ 
turbation  theory  yields  the  tunneling  time  of 
1.10x10  1Ls.  Since  the  tunneling  time  at 

Nt  =  lO^^cm  ^  is  much  less  tnan  the  boron  BE 
In 

decay  time  of  10  ®s  (Auger  lifetime)1^'17,  BE 
tunneling  will  occur  before  decay. 

We  summarize  the  results  of  the  present  calculation  by 


pointing  out  that  PL  can  be  used  as  a  quantitative  tool  to  measure 

the  impurity  concentrations  in  a  doubly  doped  semiconductor  if 

14  -3 

their  concentrations  are  less  than  10  cm  because  the  theory 


predicts  no  BE  tunneling  for  all  the  systems  considered  here. 

As  a  result,  PL  can  be  an  excellent  quantitative  tool  for  eval¬ 


uating  hiun  purity  silicon  production  and  processing  st. 


ei)S  such 


as  silane  and  polycrystalline  Si  purity  analysis,  as  well  as 
monitoring  the  compositions  of  epitaxial  silicon  films.  The 
silicon  materials  industry  is  currently  exploiting  the  PL  tech¬ 
nique  for  these  purposes. 

We  conclude  this  section  with  the  following  recommendations 
and  suggestions: 

1.  In  order  to  confirm  the  predicted  reverse  BE  tunneling 

from  a  deeper  impurity  to  a  shallower  impurity  when  there  is  a 

small  difference  in  their  ionization  energies,  we  recommend,  for 

example,  measuring  the  relative  PL  intensity  I  /I  at  ND  = 

B  Gel  B 

13—3  17  i o 

10  cm  and  N_  =  10  -10  cm  for  the  system  Si:(B,Ga)  or  the 

corresponding  Si:(B,Al)  system. 

2.  Since,  as  indicated  in  Eq.  (23),  the  tunneling  factor 
R(l,2)  strongly  depends  on  the  FE  capture  cross-section  ratio, 

1 / o 2 ,  it  would  be  very  useful  to  measure  and/or  calculate  that 
ratio  accurately,  perhaps  using  a  model  analogous  to  positronium 
capture  by  a  neutral  atom. 

3.  Historically,  minority  lifetime  killer  dopants,  e.g.  zir- 

1  8 

conium  have  been  added  to  Si  in  an  attempt  to  improve  detector 
performance  in  radiation  environments.  It  is  possible  that,  for 
excitons  in  Si,  a  deep  acceptor  could  play  an  analogous  role  in 
a  more  benign  fashion. 

4.  Since  current  silicon  growth  and  processing  technology 
has  not  yet  reached  the  point  where  all  defects  save  the  major 
dopants  are  eliminated,  it  would  be  a  valuable  extension  of  this 
work  to  consider  the  exciton  tunneling  problem  in  a  multi-component 
material.  Also,  due  to  the  increased  attention  devoted  to  Si 
counter -doped  extrinsic  detectors,  it  would  be  useful  to  evaluate 

24 


compensation  effects  by  appropriately  including  donor  bound 
exciton  capture  and  transfer  in  the  rate  equations  used  in  this 
model.  These  modifications,  while  useful  in  assessing  the 
viability  of  PL  in  quantitative  determinations  of  physically 
realistic  materials,  are  outside  the  scope  of  the  present 
investigation . 
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NOTE  TO  ACCOMPANY  FIGURES  3  THROUGH  32 


The  range  of  the  concentrations,  N^,  of  each  of  two  impuri 
ties  in  silicon  is  1012  -  1020cnf3.  The  letter  "a"  immediately 
following  a  figure  numeral  refers  to  the  two  dimensional  plot 
of  the  tunneling  factors  R(i,j)  for  a  given  free  exciton  capture 
cross  section  ratio  for  the  system  Si:(i,j)  as  a  function  of  the 
concentration  N ^  of  impurity  j  for  the  given  concentrations 
of  impurity  i.  The  letters  "b"  and  "c"  immediately  following  a 
figure  numeral  refer  to  two  different  views  of  the  three-dimensional 
plot  of  the  tunneling  factor  R(i,j)  for  a  given  free  exciton 
capture  cross-section  ratio  as  a  function  of  NT  and  NT. 


Figure  3a.  R(B,In)  for  .1  /■:  =  0.1  as  functions  of  N  for  given 


Figure  5a.  R(Al,In)  for  o  /j  =  0.05  as  functions  of  N  for  given 


R(Al.In)  for  : _ , /  =  0 . 2  as  functions  of  N_  for  qiven 
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Figure  13a.  R(Ga,In)  for  a_  /a  =  5  as  functions  of  N*  for  giver 
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l^ure  14a.  R,Ga,In)  for  '  /  -  20  as  functions  of  N  for  give 
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Figure  16c.  Another  view  of  R(B,A2)  for  j/o  =  0.02  as  a  function  of  and  N 
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Figure  21a.  R (B , Ga)  for  oB/o  =  0.005  as  functions  of  N  for  given 
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Figure  23a.  R(B,Ca)  for  Ug/o  -  1/14  as  functions  of  N.,  for  given 


Figure  24a.  R(B,Ga)  for  o/o  =  0.5  as  functions  of  N  for  given 
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Figure  25a.  R(B,Ga)  for  Jg/o,  =  2  as  functions  of  N  for  given 
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Figure  27a.  R(Al,Ga)  for  o  ,/u  =  0.05  as  functions  of  N  for  given 
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igure  27b.  One  view  of  R(Al,Ga)  for  cAl/^Ga  =  0.05  as  a  function  of  and  N 


Figure  27c.  Another  view  of  R(AJ,Ga)  for  a  ,/o  ~  *^>05  as  a  function  of  and 


ven 


Figure  28b.  One  view  of  R(Al,Ga)  for  U-./c  =  0*2  as  a  function  of 


Figure  29a.  R{Al,Ga)  for  o^,/ar  =  1  as  functions  of  N  for  given 


as  a  function  of  N  and  N 


Figure  30a.  R(Al,Ga)  for  aW  ~ n.a  =  2  as  functions  of  N  for  given 


Figure  31a.  R(Al,Ga)  for  aAj_/JQa  =  10  as  functions  of  N  for  giver* 
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